Potential field migration for rapid imaging of gravity gradiometry data
نویسندگان
چکیده
The geological interpretation of gravity gradiometry data is a very challenging problem. While maps of different gravity gradients may be correlated with geological structures present, maps alone cannot quantify 3D density distributions related to geological structures. 3D inversion represents the only practical tool for the quantitative interpretation of gravity gradiometry data. However, 3D inversion is a complicated and time-consuming procedure that is very dependent on the a priori model and constraints used. To overcome these difficulties for the initial stages of an interpretation workflow, we introduce the concept of potential field migration, and demonstrate how it can be applied for rapid 3D imaging of entire gravity gradiometry surveys. This method is based on a direct integral transformation of the observed gravity gradients into a subsurface density distribution that can be used for interpretation or as an a priori model for subsequent 3D regularized inversion. For large-scale surveys, we show how migration runs on the order of minutes compared to hours for 3D regularized inversion. Moreover, the results obtained from potential field migration are comparable to those obtained from regularized inversion with smooth stabilizers. We present a case study for the 3D imaging of FALCON airborne gravity gradiometry data from Broken Hill, Australia. We observe good agreement between results obtained from potential field migration and those generated by 3D regularized inversion.
منابع مشابه
Potential field migration for rapid 3D imaging of entire gravity gradiometry surveys
Introduction Gravity gradiometry has come to be routinely integrated into exploration workflows, since it can provide an independent measure of the subsurface density distribution. 3D density models derived from gravity gradiometry data are used to improve velocity models used in seismic imaging of complex salt and basalt structures. The advantage of gravity gradiometry over other gravity metho...
متن کاملRapid imaging of gravity gradiometry data using 2D potential field migration
We introduce 2D potential field migration for the interpretation of gravity gradiometry data. This method is based on a direct integral transformation of the observed gravity gradients into subsurface density distributions which can be used for realtime imaging or as an a priori model for subsequent 3D regularized inversion. Potential field migration is very fast and stable with respect to nois...
متن کامل3-D Cauchy-type integrals for terrain correction of gravity and gravity gradiometry data
S U M M A R Y Fundamental to complex analysis is the Cauchy integral theorem, and the derivation of Cauchytype integrals. For over 40 yr, Cauchy-type integrals have been used to describe analytical continuation, establish the location of singular points, and study non-single-valued solutions of inverse problems in 2-D potential field theory. In this paper, we revive this interesting and fundame...
متن کاملRapid gravity and gravity gradiometry terrain correction via adaptive quadtree mesh discretization
We present a method for modeling the terrain response in gravity and gravity gradiometry surveys utilizing an adaptive quadtree mesh discretization. The data-dependent method is tailored to provide rapid and accurate spatial terrain corrections for ground or draped airborne surveys. The surface used in the modeling of terrain effect at each datum is discretized automatically to the largest cell...
متن کاملGICUDA: A parallel program for 3D correlation imaging of large scale gravity and gravity gradiometry data on graphics processing units with CUDA
The 3D correlation imaging for gravity and gravity gradiometry data provides a rapid approach to the equivalent estimation of objective bodies with different density contrasts in the subsurface. The subsurface is divided into a 3D regular grid, and then a cross correlation between the observed data and the theoretical gravity anomaly due to a point mass source is calculated at each grid node. T...
متن کامل